home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / SPLINE.DEM < prev    next >
Text File  |  1991-04-29  |  769b  |  34 lines
  1. PROGRAM d3r3 (input,output);
  2. (* driver for routine SPLINE *)
  3. CONST
  4.    n=20;
  5.    pi=3.1415926;
  6. TYPE
  7.    glnarray = ARRAY [1..n] OF real;
  8. VAR
  9.    i : integer;
  10.    yp1,ypn : real;
  11.    x,y,y2 : glnarray;
  12.  
  13. (*$I MODFILE.PAS *)
  14. (*$I SPLINE.PAS *)
  15.  
  16. BEGIN
  17.    writeln ('second-derivatives for sin(x) from 0 to pi');
  18. (* generate array for interpolation *)
  19.    FOR i := 1 to 20 DO BEGIN
  20.       x[i] := i*pi/n;
  21.       y[i] := sin(x[i])
  22.    END;
  23. (* calculate 2nd derivative with spline *)
  24.    yp1 := cos(x[1]);
  25.    ypn := cos(x[n]);
  26.    spline(x,y,n,yp1,ypn,y2);
  27. (* test result *)
  28.    writeln ('spline':23,'actual':16);
  29.    writeln ('number':11,'2nd deriv':14,'2nd deriv':16);
  30.    FOR i := 1 to n DO BEGIN
  31.       writeln (i:8,y2[i]:16:6,-sin(x[i]):16:6)
  32.    END
  33. END.
  34.